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ABSTRACT 


iimeracksce, Chiyristor rectifiers are wsed as high 
current-low voltage power supplies. Under certain con- 
Gitions such rectifiers introduce subharmonic resonances. 

In this paper, a limit cycling voltage regulator 
system is simulated in the digital computer. The ref- 
erence voltage and load are varied and the system gain 
is changed for better understanding of .their effects on 
the subharmonics. A Fourier analysis of the output wave 
Shape then determined the existence and amplitude of 


subharmonic components. 
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Pot ee rOENG SYSTEMS, THEORY OF OPERATION 


A. INTRODUCTION TO LIMIT CYCLING SYSTEM 

iimonachrCe,eiG mightp be required to have a control 
system which will generate a limit cycle. In this thesis 
the goal is to investigate the characteristics of a limit 
Cycling voltage regulator. This kind of voltage regulator 
System is potentially useful as a low cost source of well 
meoulatead direct current. 

The system introduced in this paper has linear and 
nonlinear elements. One of the tools for linearizing non- 
imear Systems is the describing function method. In deal- 
mieewith Limit cycling systems the ordinary describing 
function method predicts Bie lint cycle. but 15 not ade— 
quate for a complete analysis. The Dual input describing 
iMmeculen 1s 2 useful method for investigation of the limit 


cycling system when subharmonic resonances are anticipated. 


Be DUAL INPUT DESCRIBING FUNCTION (DIDF) 

West (4), introduced a new mathematical device which 
wecaltedsihe “Dual input Describing function" (DIDF). It 
is an extension of the conventional describing function 
which linearizes the nonlinearity when forced by two inputs. 

Micsoualeinpue Describing function is a very powerful 
Laolteresthe investigation of Limit cycling systems. It 


Mem ocmuscdeLOr Ppreagictlon of subharmonic resonance. 


Tet 





For better understanding of a limit cycling system a 
general block diagram is introduced in Figure l-la. In 
this limit cycling control system N(A,w) represents the non- 
linear element, G(jw) the linear element. The DIDF is 
based on the error input to the nonlinear element. For a 


well designed system the error X(t), is equal to 
X(t) = r(t)-C(t) 


and would be small so that any limit cycling period To, can 
be modeled as shown on Figure 1-l1b. In tre noaed the sin- 
usoid amplitude A is associated with one input and the D.C 
bias amplitude is associated with the second input. Then 


the input to the nonlinear element is 
X(t) = BtAsiné 


Now, the calculation of the DIDF for the nonlinear element 
N(a,w), is accomplished by computing the ratio of the out- 
put amplitude spectrum of the nonlinearity to its input 


Watehn gives two DIDF. 


5 we 
N6A;B) ee Je 
A 
and 
= By 
N (A,B) ie 


iiemirrst nonlinearity is considered as limit cycle DIDF 
and will be used to predict the system limit cycle. The 
second DIDF is considered as a linear gain for range B/A 


and leads to linear system description due to input. The 


ike 








NONLINEAR] Y/(+ 
SYSTEM 


TTCURESEola.» Basic limit cycling system. 
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B B 
' wi 
ie 6 eran 


PEGURESI 1b. Representation of input-output relation of a 
nonlinear element. 
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FIGURE 1.1lc. Input and output wave form for an ideal relay. 
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Matmematical definition of both DIDF's follow directly from 


Fourier analysis and they can be written as 





N (A.B) = Re[N (A,B) ]+JI_ [N (A,B) ] 
27 
_ eS s 2 
man J y (B+Asiné )sinédé 
271 
bm <3 ; 
=a J y (B+tAsin6é)cosédaé 
o) (1-1) 
and 
27 
N(A,B) = + J y(BtAsind)ae 
0 (1-2) 


Equation (1-1) can be reduced since integration of the co- 
sine term of N643B) Piethtominvervals 1 9zero. 9 Then the 


Eamplified form of Nj 6AB) can be written as 
27 
N (A,B) = + f y(BtAsin®)sinode 
O om1B 
O° (3) 


In the limit cycling system B/A will be small .and an incre- 


mental DIDF can be derived as B-o 


DIDF A lim [N_(A,B)] 


rao) 
and 277 
ori Lim 1 | 
DIDF = a B20 Bas y (B+tAsiné )dé 





O 


of) 


LD 


Oe ee!” 





Cc. DIDF FOR IDEAL RELAY 
Let us consider an ideal relay with drive levels iD, 


feire L.ilc. There the input to the nonlinearity is 


X 


BtAsinw 9t 


BtAsin#é 


The ideal relay has a single valued nonlinearity, there- 
more, Computation of the limit cycle DIDF is done by using 
the result obtained in equation (1-3) and applying it to 


Figure 1.lc. Then the DIDF can be written 








271 
—_ ome e 
CB) —— i y (BtAsin6é)sinedé 
T+6y 27-6, 
= 2 J vsinedo + f (-d)sineds 
u oO T+64 
27 
+ y Dsinedé 
27-6, 
_ AD 
arn cos6, 
and 
coll 
6, = sin (B/A) 


Up 
ee. Te a ears 
N(A5B) TA 1 - (B/A)? (125) 
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The DIDF, associated with the signal is 








1 T+6y 27-0 27; 
N(A3B) aE Ddé ee Did vee 
eD 
oe 
= a sin * (B/A) (1-6) 





iieb>o, forumla (1-6) can be written 


N.(A,B) = lim [N,(A,B)] 


2D 
TA Cle, 





N.(A,B) 
Pmaamechie OID [or two sinusoids input is found as: 
X = Asin(w t+6)+Bsinw.t 


in which ae is the limit cycle frequency and Ws is the sig- 


nal frequency. Then the DIDF will be 


: 2 & 
be Vg) 1 |B 3 B 
(1-8) 
Megs) = <2 | i+ = B) ree (jie | 
a? TT 8 IA 64 {A 
(1-9) 


If this result is compared with-the approximate DIDF, it can 
Pewsoalagethay the approximate DIDF which is derived as equa- 
tion (1-5), (1-6) are valid for 5 percent accuracy but if 


the following conditions are satisfied: 


ay 





>| 
A 
LWW} 


= aoe 
A 
WI] 


O 


(1-10) 


Migem it iS possible to have 2.5 percent accuracy with re- 
Poces to the exact DIDF. 

Therefore, for the ideal relay case the limit cycle 
DIDF are always valid.to better than 5 percent accuracy re- 


ferred to the exact DIDF, under the restriction of (1-10). 


Di LIMIT CYCLING SYSTEM CHARACTERIZATION 

The system here considered is a kind of system which 
always has a limit cycle. Most limit cycling systems in- 
clude a pre-filter like that shown in Figure 1.2. The 
feneamental purpose of this filter is to reshape the command 
signal amplitude to the input such that at the input of the 
nonlinearity inequalities (1-10) are satisfied. This brings 
a design problem of the pre-filter ron will not be con- 
Saccored in this thesis. 

i Recamency Yomain Description 

To investigate sinusoidal frequency analysis of 

linmtoeecycling systems, the system model is illustrated in 
Figure 1.2. r(t) is a sinusoid and the input to the non- 


linearity N is 
X = Asin(w t+6)+Bsinw .t 


The goal here is to investigate the interaction between two 


sinusoids, through the nonlinearity and differents inputs 
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Cc tedeioa 


I 


7. 





are considered in this analysis. Figure 1.3 shows the range 
of possibilities of the signal and limit cycle. For PFig- 
ure 1.3a the output through the nonlinearity is nearly sin- 
Usoidal at the limit cycle frequency but is phase shifted 
with respect to the limit cycle. It can be said that sig- 
nal through relay is zero which means N.=0. Also it may 
be said that for signal frequency greater than the limit 
cycle frequency but satisfying the amplitude ratio con- 
dition also gives N.=0. 

ya) Figure 1.3b the amplitude ratio is satisfied 
but the frequency ratio is 0.5. In this case the signal 
through the nonlinearity would not lose identity and the 
amplitude ratio only will depend on the piasessitit of the 
input component to the nonlinearity. Then for the case 
shown in Figure 1.3b, it can be concluded that a eee 
horse Lrequency ratio equal to 0.33, N, affects the input 
signal amplitude by factor of 16 percent and a phase of 
4 degrees. For a frequency ratio greater than 0.33, (but 
amplitude ratio satisfied) Nye can be Monerriened independent 
of relative time position since the change in amplitude and 
phase ratio are never greater than 1 percent. Therefore as 
eeresult of these Conclusions the open loop transfer function 
can be written as 


Ww 
& . =O. 
GEG JOIAN, -Hg-@(Je) for |w|< 3 


Ww 
O 
jw)= = < < 
G4, (5) YN. -Hog -G(jw) for 3 <|w|<w, 


Go, (5w)=0 for w 9 < lo (1-11) 
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PGUiREeie. sa. AR = 0.33, FR = 1.00 


%= a+b 
a 
b 
t 
PACURPIesben AR = 0.33, FR = 0.5 
x=a+b 
ni 





PECURE ese eR = 0.55, FR = 0.16 


FIGURE 1.3. Range of possibilities of input. 
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where y represents a complex number depending on the rela- 
tive phase of these signals. 
2. Pole-Zero System Characterization 
Frequency domain analysis shows that in the range 


between -W/3<w sss rad/sec along the imaginary axis, if 


Ss 
the amplitude ratios are also satisfied the system behaves 


linearly. In the ranges w /3<w <0, and -W Sw <w/3 rad/sec 


S S 


Dew sienal gain of the nonlinearity is a function of time, 
Pmcenor the remaining portion of the jw axis, the signal 
gain is always zero. Now, consider Seenent tals for AE 
the amplitude ratio is satisfied, (amplitude ratio refers 

to ratio of peak exponential to peak limit cycle amplitude) 
it may be noted that the exponential excited output of a 
relay is relatively independent of the time origin of the 
input exponential signal in the existence of the limit cy- 
liermenemcnewstenal lasts £or two or more limit cycle periods 
Figure 1.4 shows this situation. Then, the minimum value 


Si svhne input time constant is, 
Tt Tae ~ 2 | 
min W) 


1 0. 
ee S 


min 





Then, the pole-zero configuration of the system may be 
Gescribed. Consider the complex plane "p" to be divided 
into three regions, shown in Figure 1.5. 

The nonlinear element N of the system will be 


considered as a gain element and N is derived for the relay, 


ee 





a+b 


TOTAL RELAY OUTPUT 





LIMIT CYCLE PORTION 





EXPONENTIAL PORTION 





FIGURE 1.4. Response of ideal relay in the presence of 
Dimi cy Cile:, 


ZS 





ae 


_— Jw, 
7 X% OPEN LOOP POLE 
CL} CloSED LOOP POLE 


PLGURES 5. Root locus plot for G(p). 


24 





PreGetore, it Will not contain any phase shifting. If the 
open loop system gain is adjusted to one-half of the value 
@eevne limit cycle fain, it is possible to automatically ob- 
maim a limit cycling system. The closed loop roots of the 
System can be approximately determined, such that, if the 
gain at the crossover of the imaginary axis NEN which is 
eae limit cycle gain, and at gain value of Notes = one 
of the closed loop roots will be close to the intercept 
point of the w /3 eirelle on the real axis which is shown as 
PA I a and the other complex pair of roots would be 
very close to the circumference of the Wo circle. This 
result determines the dynamic response of the limit cycling 
Syoucm. 

One Of the important points about limit cycling 
Syeuems is from the stability point of view: if the sensi- 
tivity of the system is increased then the limit cycle am- 
plitude increases but the gain Ne decreases and always the 
S.N. DeOGdilc oe remaincweonstanyv. Closed loop root locations 
do not change, only the limit cycle amplitude will be 
changed. Also, it is possible to have instability in the 


system but may be the system response does not change. 


Ee CEOSED TILOOP PROPERTY 
Consider a limit cycling system as shown in Figure 1.2. 


Then the open loop transfer function of the system is 


i 


Cre 8), 


aa N-H,(p)-H5(p) .G(p) 


NoHo GCP) Se: 


2S 





The DIDF for the ideal relay were found in equations (1-5) 
and (1-6) and they may be approximated as, 
~ AD 
No * TA 
2D 
s TA 
emaethe ratio of these two, 


N 


Jes 


N 
ie 


hole 


then N. may be written in equation (1-12) 


Gqi Jw) Ne Hog Gliu) 


aq Ne tecde 0”? 


Mamorder to sustain the limit cycle 


NHR GC5e) =- i pedis) 
muse be satisfied. Then, 
fiacs ae 
Seg I = esti) 


tess a kind Of proof of the explanation given earlier 


heomeeopen 1600 gain of limit cycling systems. 
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ies PV OnMANCE CHARACTERISTICS OF THE LIMIT 
CYCLING VOLTAGE REGULATOR 
A. SYSTEM DESCRIPTION 

The block diagram of the voltage regulator system used 
miechiis thesis is shown in Figure 2.1. The complete circuit 
diagram is shown in Appendix A. 

Ihe firing circuit, instead of using a separate cir- 
euice for each thyristor uses a relay to turn all thyristors 
on or off. If the relay is on a positive voltage is applied 
morali- thyristor gates and this causes GCOnduction in the 
maettiter seit the relay is off then the rectifiers wiil 
mor conduct. . 

The rectifier unit is a three phase full-wave recti- 
mie~emade Of a combination of six thyristors. The thyristor 
PeeoesOlMidnsvuabe Geyice which is e€quivalent to a thyratron. 
One thyristor device is a silicon controlled rectifier. It 
is a four layer device and may be likened to a PNP transis- 
more joined to an NPN transistor and sharing:a common col- 
lector-base junction. The thyristor and the characteristic 
curve of the thyristor are shown in Figure 2.2. 

As a load unit, a second order filter with very low 
Cemipaient aeter 1s Used. It is an LCR filter. The values 
eel. c are Consvent and variation of the damping ratio is 
obtained by varying the load resistance. 

The system has a nonlinearity due to the relay and 


the thyristor characteristics provide a second nonlinearity. 


27 





AC SUPPLY 


ee THYRISTORS FILTER 


CIRCUIT 





FIGURE 2.1. Basic block diagram of the limit cycling 
voltage regulator. 
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These two nonlinearities are treated as one nonlinearity in 
analysis of the system, providing the simplified block dia- 


gram of the system as shown in Figure 2.3. 


B. STEADY STATE ACCURACY 

It is important to investigate the system response to 
aqiiferent inputs. For this purpose the voltage regulator 
system is analyzed for a load resistance R=0.82, and dif- 
ferent reference voltages applied to the system. The im- 
portant feature to be studied is the output characteristic 
of the system. For different reference voltages the out- 
puts obtained were subjected to Fourier analysis. Table 2--1 
Shows the result of this analysis. in Table 2-1 results 
taken from digital computer calculation of the Fourier se- 
ries and subroutine FORIT used for this analysis. In the 
Table 2-1 A ie the DC velue of the output. 

Mhese results show that the system output is not equal 
to the required value. Then the difference may be defined 


as a, 


For high reference voltages the system has large a and if 


ef 


eomcianeteseit' s sign, that is, for low voltage cases 


a, gradually decreases, a also decreases. At Me = 52 


MicmovovemmouupULeVvOlbaze 15 4a little bit higher than Leyak 


This may be good but when us decreases the system ripple 


ef 


instability increases which is not desired for the system. 


The numerical results show that if Ve is decreased gradually, 


ef 
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BEGURES ZS oilmplicaed block diagram of the limit cycling 
voltage regulator. 
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FOURIER ANALYSIS FOR R=0.82 


REFERENCE LIMIT CYCLE] SUBHARMONIC] LIMIT CYCLE | SUBHARMONIC] ORDER oF 
VOL TAGE A, (vor) AMPLITUDE | AMPLITUDE | FREQUENCY | FREQUENCY | suBHARKONIC 
(Hz) 


eA 45.0 Wz 


0. 4 
48 Aes ted. 505) 0.475 56.2 281 Ge 
ig . ue 
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subharmonics become more effective and this may make the sys- 
tem unstable. For Vio» = 58 volt, the amplitude of the sub- 
mermonic is quite small (about 0.146 volt) but for ae. = 

46 volt, the subharmonic amplitude becomes 0.635 volt which 
is quite large compared with 0.146 volt. Also the results 
show that the limit cycle amplitude shows an opposite re- 
action. Increase in the subharmonics amplitudes caused a 
decrease in limit cycle amplitude. Figure 2.4 shows the 


output of the voltage regulator at Vip = 60 volt and 


Weer = 58 volt. 


C.. TRANSIENT RESPONSE OF THE SYSTEM 

DOimnvesta*avdon Of the transient response of the sys- 
tem only the linear parts are taken into account. The be- 
Maer eOl che Nonlinear part is not considered. The linear 
jmieer Ol tne System is assumed to be as shown in Figure 2.5. 


And the transfer function of the system is 


G(s) = Vout _ Te/AL6 
Vin S- FIR OSU Le) (=) 





In general second order system characteristic equations can 


be written in the form 
S*+2&w Stw 7= 0 
n n 


then for the filter, 
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PiGukimzesmebitcer and load of the system. 
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Purine analysis of the system L and C are held constant and 


only the value of R is changed 


C= On it, 
L = 0.0005h. 
therefore, 
One Tie rad/sec 
emoetor R = 0.82 
2 = 0), OL aieyn 


which is a very low damping factor. 


Miewseuvtling time of the filter is 


vous oe 
n 
ier, 
Dee = 2.5 sec for R = 3.22 
t. = 1.23 sec for R = 1.62 
te =sOnous5esec for Rk = 0.62 


The filter transfer function can be written in Bode form, 


, a ill 
G(Jw) = LC( jw) *+L/R(jw)+1 


Madetor RnR = O.o82 


ne) co 
el eee Gi 2 Ga) (2372 10" (ju) 


(2=2} 
The phase and amplitude curves of the Bode diagram are 


drawn with the IBM-360 digital computer and the curves are 
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illustrated in Figure 2.6. The analysis of the closed loop 
Sem@earacteristic of the filter is shown on the Nichols plot 


iorigure 2.7. 


»'. PREDICTION OF LIMIT CYCLES 

Percepredrccion of the limit cycle frequency, the basic 
@escribing function method is used. The describing function 
for the nonlinearity of the system which includes relay and 
rectifier is derived by Leszcynski (Ref. 15). And for the 
linear part of the system for this prediction equation .(2-2) 
is used. Figure 2.8 and Figure 2.9 shows the Nichols plot 
of the limit cycling voltage regulator for ee 64 volt 


and We —oOcevotl bothoat RR = 0.60. Figure 2.8 gives the 


ef 


Mimat cycle frequency Poo =—ss5i24 aad Jimit cycle amplitude 
B= 3.6544 volt. This agrees with the experimental result 
shown in Table 2-1. Also for Vp = 62 volt the Nichols 
chart of Figure 2.9 gave the limit cycle frequency Pi StH2 
and limit cycle amplitude about B = 2.4632 volt. Experi- 
mental results in Table 2-1 gives f,, = 36.4HZ and limit 
eycle amplitude about B = 2.82 volt which agree fairly well 


with the result obtained by Nichols plot. 


De SUBHARMONICS 

One of the purposes of this paper is also to investi- 
gate subharmonic resonances in the system. Table 2-1 for 
Fourier analysis shows relationships between Maer and the 
subharmonic amplitude and frequency. Results show that 
changes in eee affect subharmonic amplitudes. If ee is 


increased the subharmonic amplitudes also increase and for 
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FUGUIRESe2 0), 
diagram of the 
Hincar part. 
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PmecuRE 2.7. Nichols plot of the filter for R = 0.82. 


59 


&)= 120 





G(jw) 


udeldl.S wde=l40 ; 


20 U) 2138 


a) = 145 


—-* 
W 


g wW=153 


10 


wn 


(DB) 


-10 





+150" -180 -150° -120 SN 


FIGURE 2.8. Nichols plot of system for ees GANG R=0 On? £16739: 


40 





20 


oO W=186.5 


(DB) 





° 


+150° -180° -150° -1l20 ~90 


PUGUR EZ. ose Nichols plot of the system for V =62V., R=0.82, 


fie = I alvin 


ref 


41 





V < 50 volt subharmonics become very important. System 


ref 
mehavior £ives a ripple instability at the OuvpDUuL. this n= 
Peorlity May De Caused by the high amplitude subharmonics. 
Pueure 2.10 shows this kind of output for aes = 48 and 


Uo 46 volt. 

From the frequency response point of view, it may be 
eencluded that changes in ar CG, aitect sappreciablly the 
frequency of the subharmonics. The basic frequency of the 
Subharmonic tends to remain at 25 cps. Table 2-1 shows that 
changes in the system load Be rance: R = 0.82, R = 1.62, 
R = 3.22 caused no change in the subharmonic frequency and 
io like manner changes in ume reference voltage did not 
Seanee the subharmonic Trequency. These results shown in 
Table 5-1, Table 5-2, and Table 5-3. It may be concluded 
tiao any change in the system will not effect the subhar- 
monic frequency since the system forced the subharmonics 
toward the 25 HZ value. But the amplitude of the sub- 


harmonic is quite sensitive to any change in the system or 


any disturbance to the system. 
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TIil. INTRODUCTION TO THEORY OF SUBHARMONIC RESONSNCE 


Subharmonic resonance is a phenomenon which occurs in 
nonlinear systems. Usually the nonlinear system subharmonic 
resonance frequency is a fraction of the input frequency. 
As an example consider the limit cycling voltage regulator 
Meeeme which 1s introduced in this thesis, or a relaxation 
oscillator, or a moving-coil loudspeaker with a non-uniform 
magnetic field are the common examples of the subharmonic 
resonance. This phenomenon can also occur in servo mecha- 
te Mmomvleiesaluration Or With other nonlinearities. In all 
of these cases the lowest frequency present at the output 
Pwo oupimliil tiple of the input frequency and usually aor 
Mearetne Mavural irequency of the system. The system out- 
pureamoligude would be different depending on the existence 
PeeeneesUbmaarmonic. Such Subharmonics can increase or de- 
crease the output amplitude. For systems in which the out- 
put is intended e reproduce the input as closely as 
Pessioleyeaseis the motion of Che loudspeaker cone or a 
servo motor, this phenomenon can introduce an intolerable 
Ceisvorey lon nh lSO 4b 1S observed that in limit cycling vol- 
tage regulator systems, subharmonics can cause instability 
Otetune system. Also one of the important points, subhar- 
Mmonmire resonance does not occur in ae linear systems. If 
subharmonics are not wanted the best way to eliminate them 
is by changing the parameters Ousuhewinear sysvem or plac- 


ing a compensation element rather than trying to change the 


nonlinear part of the system. 


yy 





iGei1s possible to predict subharmonic resonance in 
the nonlinear Soto yeUcaneg tie Dual input Describing 
function and by solution of the nonlinear differential equa- 
ti20n. There are some disadvantages to the second method. 
Since the output of the nonlinear element is expressed as a 
power series of the input, successive approximations to the 
Pemeecs SOlUtION converge rapidly only if the amount of 
nonlinearity is small. And also the solution of the dif- 
Beeenuial equation does not give a physical picture of the- 
mechanism of subharmonic resonance. 

Subharmonic resonance normally exists only in lightly 
damped systems with nonlinear energy delays or with nonlin- 
ear restoring forces. In systems giving rise to subhar- 
Homies. cne Nonlinear element is frequently characterized 
by hysteresis. Typical nonlinearities are backlash, mag- 
Hewtewhysoteresis loops and relays, etc. It is also impor- 
Wale to) note that generation of subharmonics is dependent 
Mimo ueraleCOndiutons, 1,.€6€., they may occur in a certain 
range of amplitudes and frequencies. If ina system a 
subharmonic is obtained it stays “locked in" over a reason- 
Somicmrregquency, range. If the driving frequency is changed 
to a new value, the system output frequency changes also 
to a new value which is w/n, w being the driving frequency 


and n order of subharmonic. 
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IV. DUAL INPUT DESCRIBING FUNCTION FOR 
PREDICTION OF SUBHARMONICS 

In the voltage regulator system, power supply ripples 
feed around the loop and interact with the nonlinearity of 
Wiemsysuem and cause instability of the system. The limit 
eucle resulting from this instability is at a submultiple 
of the ripple frequency and is commonly called a subhar- 
monic resonance. Normal describing function theory cannot 
Peeguver vile Existence of the er neunontc resonance. But 
temusiane the Dual input describing function method, i is 
possible to find a mathematical description of the nonlin- 
Cam@epart Gf the system which can be used to predict the 
existence of these resonances. In general the forward 
Peaniecter tunecvion of the voltage regulator system can be 


written 


Vout = N(jw) .@(jw) 


Cad); 


and 


N(jw) = No (jw)+N, (jw) 

ioewhaech N (jw) is the DIDF due to the DC input to the non- 
linearity of the system and N, (jw) TomauetoenG input Co 
the nonlinearity. G(jw) is the transfer function of the 
linear part of the system. If E represents an input to 

the nonlinearity then equation (4-1) can be rewritten, 


Vout _ Nj GCJw)tn, (jo) 
: Gee) 
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Frequency analysis is accomplished in the usual fash- 
muwmmeoy applying the Nyquist criterion. The result could 
be interpreted as an encirclement of the -ltjo point. 
Therefore, equation (4-2) can be drawn in vector form as 
Seeowo in Pigure 4.1. It is seen from Figure 4.1 that the 
N,-G(jw) vector is placed at the tip of the NGC jw) vec- 
tor. Because N (5) meopiresents the DC part of the DIDF 
and Nj (J) represents the AC part of the DIDF, N (5) will 
have one component but Nj (jw) will have many components as 
the phase angle changes. But during this variation w = OS 
must be constant. Then if (for different w) this vector 
plot is repeated Figure 4.2 would be obtained. The N,-G(jw) 
Mmeerer would be placed on the tip of the NjG(5o) as phase 
iemeeciorcd,wepy LOollowing this procedure many curves can be 
drawn and if any of these curves N,-G(jw) encircles the 
-l+jo point then in Nyquist sense the system may be un- 
stable. Such a system may have a subharmonic resonance, 


but this is not guaranteed. 


an 








(N+ N,). G(joy) 


FIGURE 4.1. Vector diagram using DIDF. 
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PIGURE@4=2s5 Analysis of the DIDF. 
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eon omOn tne se XPERIMENTAL RESULTS 


A. ANALYSIS OF LIMIT CYCLING VOLTAGE REGULATOR 

line @ealeas ve find out the exact cause of subharmonics 
fimeode limit cycling voltage regulator system. For this pur- 
pose the IBM-360 computer was used for solving system dif- 
ferential equations and subroutine INTEG 2 is used. Appen- 
dix B shows this computer program and as an example one of 
the outputs. The output graph was obtained for three daif- 


ferent loads and in each case Ve is changed from 64 volt 


ei 
to 46 volt. Each run was permitted enough time so that the 
system reached steady state. From these output plots by 
analyzing the wave shape of the output, the results in 
Table 5-1, Table 5-2, and Table 5-3 were obtained. 

Those tables showed that the system has mostly 1/2 
order subharmonics and a few 1/4, 1/8 subharmonics. Then 
iG may be said that the limit cycling voltage regulator 
system has an even order nonlinearity. Because most ex- 
perimental results have shown that nonlinearities of even 
order give even order subharmonics. 

In the Tables amplitudes are given peak to peak in 
volts. These amplitudes are given as the observed values 
of the output's highest spike and lowest spike. These low- 
est and highest spikes become Barer ti tae load resistance 


issimerecased. Imcreasing load resistance caused increase 


in the ripple instability at the output. For many cases it 
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ANALYSIS OF OUTPUTS FOR R=0.8 2 





TABLE 5-I 


sul 





ANALYSIS OF OUPUTS FOR R=Il.62 


REFERENCE LARGER SMALLER |LIMIT CYCLESUBHARMONIC 
VOLTAGE | AMPLITUDE [AMPLITUDE | FREQUENCY] FREQUENCY ISUBHARMONIG 
Vaer (vout) 





* Subharmonic is determined by Fourier analysis. 


TABLE 5-2 


BZ 





ANALYSIS OF OUTPUTS FOR R=3.2 % 





SSubmarmonie as determined by Fourier analysis. 


TABLE 5-3 
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Was not possible to determine the order of the subharmonics 
by inspection of the output curves. In these cases a 
Fourier analysis was made to establish the subharmonic fre- 
quencies and amplitudes. 

Figure 5.la to 5.1d shows the wave form at the filter 
element and from Figure 5.2 to Figure 5.34 shows the output 
Pemipuceplou CO tne filter. If one of the input plots is 
Carefully observed one notices the difference in width of 
Eaese pulses. The system has sinusoidal input to the thy-~— 
mesuors and ne OucouL of this Pee ie goes to the fil- 
Pweioicewill charge the capacitor at the filter and when 


the output reaches ve the relay would be open but at this 


er 
moment the rectifier will still continue to conduct and 
Guring these times the capacitor will charge. Afterwards 
1G will discharge till the output voltage reaches Mee 

If the output voltage is much higher than ee then the ca- 
Pecmvor will adiseharge for a longer time, therefore, the 
system would give a different pulse width due to this ef- 
fect. This can be seen by the different number of ripples 
at the top of the input voltage wave form. 

Ciepucmplolvo trom Bigure 5.2 to Figure 5.11 are for 
load resistance R = 3.22. At lower voltages deciding the 
order of the subharmonic by inspection becomes impossible. 
Those outputs were analyzed in Fourier point of mien and 
those outputs obtained by using Fourier analysis are desig- 


nated with asterisks in the Table. 
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Peeuree 5.12 to Figure 5.14 shows the input plot to 
the filter with average value of the pulse width approxi- 
mately corresponding to the limit cycle frequency. 

Output and input plot for the R = 1.6 case are 
pEmowo in Figure 5.15 to Figure 5.26. For R = 0.88 case 


they are shown in Figure 5.27 to Figure 5.34. 


B. ANALYSIS OF THE D-C SYSTEM 

In this paper, one of the important things to be deter- 
mined (if possible) is the exact cause of the subharmonics. 
Therefore, to get results it is necessary to make some 
change in voltage regulator system. Then it will be easy 
momeceythe effect of the element, (i.e., relay, thyristor). 
for this purpose the rectifier part was replaced with a 
D-C source, and the system was as shown in Figure 5.35. 

For analysis of the system the C, -5000 analog com- 
puter was used and the system output obtained for R = Ore Ge 
helo, and R = 3.28%. In this system the filter unit 
Meocmene same as Tor all other tests, the only change was 
replacing the rectifier with a 65 volt battery. Fig- 
ure 5.36 shows the analog diagram for the C,5000. Ourpuc 
and input plots for different loads are shown in Figure 
5.37 to 5.39. These outputs show that the system did not 
have subharmonics. Also these runs were made for different 
Noe conditions but the system never produced subharmonics. 
This means that if the system has a nonlinear element but 
does not have a periodic forced input then the system 


never gives subharmonics at the output. In order to have 
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FIGURE 5.26. Voltage into filter, V,,, = 60V., R = 1.69 
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FIGURE 5.353. Voltage into filter, Nee = 60V., R = 0.82 
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subharmonics, a control system must satisfy two conditions, 
ive., the system must have a periodic input forcing func- 
tion and a nonlinear element. From this it might be said 
that the exact cause of subharmonic in the limit cycling 
voltage regulator is due to the nonlinear part together 
With relay and rectifier and also the periodic input to the 
system from the AC ripple. It is important to remember 
That in order to have subharmonics the system usually must 
be lightly damped. 

The block diagram shown for D-C case at rabies 5G Ss 
detail of the elements are explained in Appendix C. The 
Svouco and Input plot of system shown from Figure 5.37 to 
Figure 5.39 are for steady state OTE EAGT Figure 5.40 
shows the change in pulse duration of the relay upon the 


system reaching the steady state. 


Ge GAIN CONSIDERATION 

Momomilasceamvyestication of the limit cycling voltage 
regulator, the gain was changed from 1 to 10,000 and the 
output plot was obtained for different values of gain. 
Figure 5.41, Figure 5.42, and Figure 5.43 show these out- 
put plots. Output shows that change in system gain did 
Motsettcets bie Output wave form and it did not effect the 


subharmonic resonance. 
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PIGURE 5S.4]1a. Output plot, Vee = 60V., R = 3.22, Gain = 3.0 





5 
0.099461 0.12 Oe : ‘ a2 


FIGURE 4.41b. Output plot, Vj, = 60V., R = 1.62, Gain = 5.0 
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FIGURE 5.4lc. Output plot, Wee = GOVemeRe=) 0.00, Gain = 5.0 
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PUGURE S.42Za. Output plot, Vee ="OOV., R = 3.28, Gain = 1000-0 
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FIGURE 5.42b. Output plot, V,,, = 60V, R = 1.62, Gain = 1000.0 
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PiGURE@oe+ze. Output plot, Mee = 60V, R = 0.82, Gain = 1000.0 
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PRGURE 5.43a. Output plot, ae =SoUV,. KR — Sze, Gain = 105000 
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FIGURE 5.43b. Output plot, Vj, = 60V, R = 1.6%, Gain = 10,000 
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PUGURE S.45c. Output plot, Mae = 60V, R = 0.82, Gain = 10,000 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Unaderesonemvoda conditions a limit cycling voltage 
regulator system gives subharmonic nec onences at the out- 
Poo. this is not permissible in many applications. The 
goal of this thesis is to investigate the voltage regula- 
foeresystem and to try to find out the causes of the sub- 
mermonie resonance. 

Removal of the rectifier and substitution of a bat- 
tery removed two basic characteristics from the system: 

1. The ripple at frequency 360Hz due to the recti- 
miecation of the AC supply. 

2. Teas eed shut-off (i.e., continued conduction 
after the error detector signal gave the "off" signal) 
auewue Lhe thyristor firing characteristics. 

Since the subharmonics are observed to be submulti- 
ples of the limit cycle frequency, therefore, item (1) 
cannot be the cause of such subharmonics. Item (2) 
Seccurceab the limit cycle frequency, and results showed 
that the input voltage to the filter had a variable pulse 
duration for rectified AC but there was no pulse width 
modulation which probably caused the subharmonics. It 
Ganmoemseen trom the input plots that the ripple begins 
at different voltage levels even though the system load 


resistance and reference voltage are held constant during 
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one analysis. The thyristor characteristics apparently 
Cause this pulse width modulation. 

Another important thing is the prediction of sub- 
harmonic resonance by using the eat ipuleg er Describing 
function. If the Dual input Describing function is derived 
then by designing a compensator the subharmonics may be 
removed from the system. As a further study, the Dual 
input Describing function of the nonlinearity should be 


derived, taking into account the pulse width modulation. 
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APPENDIX B: DIGITAL COMPUTER PROGRAM FOR LIMIT CYCLING 


VOLTAGE REGULATOR 


The constants C(I) represent 


C(1) 
C(2) 
C(3) 
C(4) 
eiG>) 
C(6) 
C(7) 
C(8) 


Maximum value of input sine wave 

Ver vemoteune,  filGer inductance is henries 
Value of the filter capacitance in farads 
Veltemolloa@eresistance in ohms 
Reference voltage 

Gain of amplifier before dead zone 
Magnitude of dead zone 


Gain of amplifier after dead zone 
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APPENDIX C 


SUMMER AMPLIFIER: 
Amplifiers on the patch board are addressed as an 
even number AOOO to A126, and are represented with a gain 


symbol at the input to the amplifier 


INVERTERS : 
Sign changes are addressed from A200 to A217 and are 


represented as 


eNEUH CUREUTS 


INTEGRATOR: 
These amplifiers are addressed on the patch board as 
odd numbers (A001 to AOH7) and are integrators. They are 


represented with a symbol 


AOO7 and AO11 were used. 


LOL 





COMPARATOR: 
meerindeom device wiieh gives Ov. or -6v. logic level 


according to the comparison of the inputs. Represented 


with symbol 





Such that if 


a" 
if 


-~6v. 

(oh 

X+Y >0 < — 
C = Ov. 

(o} 
C = Ov. 

a 

X+Y<o — 
C -6v. 

a 


DIGITAL/ANALOG SWITCH (DA): 


It is a switch and is represented with a symbol 


je eocie 
SIGNAL 





The switch is closed if logic signal is -6v. and open if 
logic signal is Ov. These logic levels usually come from 
a comparator. Switching time is approximately equal to 


1 usec. 


102 





10)s 


A. 


ven 


LIST OF REFERENCES 


LivuetwomminGmorown, RNR. G&., Analysis and Design of 
Feedback Control System, McGraw-Hill, .1960. 


ijcletemcmemand Pastel, M. P., Analysis and Design 


of Nonlinear Feedback Control Systems, McGraw-Hill, 
1962. 
Gelb i eomcevemder Velde, W. E., Multiple-Input 


Describing Functions on Nonlinear System Design, 
McGraw-Hill, 1968. 


WSS@emdemcr, Analytical Techniques for Nonlinear Con- 
trol Systems, VanNostrand, 1960. 


HoveclieeerneNOnlineer Oscillation in Physical Sys- 
tems, McGraw-Hill, 1964, 


Poole mor eanarackson, FP. A., Electrical Control 
mise eri nee Vols 1. London TLIFEE, ees 


Cet [he Analysis and Design of Limit Cycling 
Ndapuivesi\ubpemante Control Systems", Doctoral 
Dassertation, MIT, 1961. 


Mavic ml ne nmmcass, A. ©., Thaler, G. J., and 
Tribes e.. “Stability Analysis of Phase Con- 


trolled Regulators", IBM Machine Technology 
Report, 14 May 1963. 
Veo mca a DOUCe. J. lb., and Livesley, R. K., 


ice lisielapUbeDeocribing Bunetion and. Its Use in 
the Analysis of Nonlinear Feedback System, Proc. 


CHa ee olen iOS Spars Beep. 4635 July 1955. 


West, J.C. and Douce, J. L., Mechanism of Sub- 
harmonic Generation in a Feedback System, Proc. 


Otminnn ey orm 1102, pe 563, July 1954. 


Cipsenned. &., and Sirdhar, A New Dual-Input 
Deoemupine Funetion and an Application to the 


Duo y sou horcea Nonlinear Systems, Trans. AILEE, 
p. 65-70, May 1963. 


Gibson, J. E., A New Dual-Input Describing Function 


aniGdwan nop lication to the Stability of Forced 
Nonlinear Systems, Trans. AIEE, p. 1-6. 


10 3 





MS 


14, 


ae 


EG, 


CoMOpe mclamvander Velde, W. E., On Limit Cycling 
Comtrpotmeoystem=s. crans. IEEE. Autom. Control, p. 
142-157, April 1963. 

Celi woe ice Oyneamnie Input-Output Analysis of Limit 


GC ClinemcConvroleoystems, Froc.”JACC, p. 3.3-l-— 
3.5-li-June nog. 
HeoZcvislaav denn atlalysiS of a Forced Limit-Cyclin 
Regulator, Thesis, Naval Postgraduate School, 1968. 


WwCc wie or ey inVveStization of Subharmonic Ripple 


Mime erecomlamne- Cycling Regulator, Thesis, Naval 
Posteraduate School, 1968. 


10 4 





ENTTIAGL DISTRIBUTION LIST 
No. Copies 


Defense Documentation Center 2 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0212 2 
Naval Postgraduate School 
Monterey, California 93940 


Wentz Kaw KS aL 
Persone Beitam sb. Mudurltigu 
Ankara, Turkey 


Deniz Harb Okulu K. Ak 
Heybeliada, Istanbul 

Turkey 

Dr. G. J. Thaler (thesis advisor) 10 


Department of Electrical Engineering 
Naval Postgraduate School 
Monterey, California 93940 


Istanbul Teknik Universitesi 1 
Blektrik Fakultesi 

Taskisla, Istanbul 

Turkey 


tices eA li GU ler 4 
Zulalecesme sok No. 2/3 
Fatih-Istanbul 

Turkey 


105 





UNCLASSIFIED 


Secunty Classification 





DOCUMENT CONTROL DATA-R&D 


(Security classification of titfe, body of abstract and indexing annotation must be entered when the overall report is classified) 


ORIGINATING ACTIVITY (Corporate author) 28. REPORT SECURITY CLASSIFICATION 


Naval Postgraduate School UNCLASSIFIED 
Monterey, California 93940 2b. GROUP 


REPORT TITLE 


Subharmonic Resonance in a Limit Cycling Voltage Regulator 


DESCRIPTIVE NOTES (Type of report and, inclusive dates) 


Master's Thess, December 1970 


. AUTHOR(S) (First name, middie initial, last name) 


Meets Guler 


REPORT DATE Ja. TOTAL NO. OF PAGES 7b. NO. OF REFS 
De eCemoe OTK 107 16 


8a. CONTRACT OR GRANT NO. 94. ORIGINATOR’'S REPORT NUMBER(S) 
- PROJECT NO. 


9b. OTHER REPORT NO(S) (Any other numbers that may bu assigned 
this report) 


. DISTRIBUTION STATEMENT 


This document has been approved for public release and sale; 
tan aiscurioutLion is unlimited. 


» SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 


Naval Postgraduate School 
Monterey, California 93940 


. ABSTRACT 


iMmpracnacessuhyristor rectifiers are used as high current- 
liemepower supplies. Under certain conditions such rectifiers 
introduce subharmonic resonances. 


ig@ethdcepe vers a limit cycling voltage regulator system is 
simulated in the digital computer. The reference voltage and 
load are varied and the system gain is changed for better 
Migerscanding ot their effects on the subharmonics. A Fourier 
analysis of the output wave shape then determined the existence 
and amplitude of subharmonic components. 


DD .2"..1473 (Pace 1) 


S/N 0101-807-6811 106 Security Classification 


UNCLASSIFIED 


A-31408 





PNCLAS SPE LED 
Security Classification 
KEY WORDS 








Subharmonic Resonance 
mma c Cycling 


woal input Describing Function 


DD J2%..1473 (sack) UNCLASSIFIED 


— J —E— 
S/N 0101-807-6821 Oe]. Security Classification A-31409 
























Thesis 


686366 
Cat 


Thesis 1256868 
86266 Giller 
c.1 Subharmonic resonance 

in a limit cycling vol- 

tage regulator. 








